Orbital polarons and ferromagnetic insulators in manganites 
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We argue that in lightly hole doped perovskite-type Mn oxides the holes (Mn 4+ sites) are surrounded 
by nearest neighbor Mn 3+ sites in which the occupied 3d orbitals have their lobes directed towards 
the central hole (Mn 4+ ) site and with spins coupled ferromagnetically to the central spin. This 
composite object, which can be viewed as a combined orbital-spin- lattice polaron, is accompanied 
by the breathing type (Mn 4+ ) and Jahn- Teller type (Mn 3+ ) local lattice distortions. We present 
calculations which indicate that for certain doping levels these orbital polarons may crystallize into 
a charge and orbitally ordered ferromagnetic insulating state. 



I. INTRODUCTION 

The interplay between the magnetic and electric prop- 
erties of Ri- x A x Mn03 (R is a trivalent rare-earth ion 
and A is a divalent alkaline-earth ion) has been ex- 
tensively studied in the light of basic physics as well 
as their technological importance.EJn La/Sr and La/Ca 
based compounds Lai-^Sr^MnOs and Lai_a;Ca x Mn03 
are metallic below the Curie temperature for 0.15 < x < 
0.5. While the transport properties of the La/Sr system 
with x — 0.3 are well-described by the conventional dou- 
ble exchange theory,BB the La/Ca system (x ~ 0.3) shows 
the colossal magnetoresistance (CMR) behavior which, 
cannot be explained by the double exchange mechanism.cl 
Neutron diffraction and x-ray absorption measurements 
of the La/ Ca system indicate that the small magnetic po- 
larons, in which spins are fecromagnetically aligned, play 
a role in the CMR effectJafl Recently, it was suggested 
that phase separation between the ferromagnetic metal 
(FM) region and the antiferromagnetic insulator (AFI) 
region is essential for the CMR behavior!] and, actually, 
the coexistence of-,the FM and AFI clusters has been ob- 
served in La/CaJ3 On the other hand, Pri-^Ca^MnOa 
is insulating for all x\ it has the AFI phase for 0.3 < 
x < 0.5 and the. ferromagnetic insulator (FI) phase for 
0.1 < x < 0.3opa although the structural details are still 
not elucidated. til It has been reported that in the Pr/Ca 
system, while the magnetic moments are almost ferro- 
magnetically aligned for x ~ 0.2, they are canted for x 
~ 0.1 and 0.3, namely, near the boundary between the 
AFI and FI phasesBEEl Recently, it has been found that 
the La/Sr system of 0.10 < x < 0.15 exhibits a FI behav- 
ior at low— temperature and may have charge and orbital 
ordering .113113 This all raises a question as to the origin 
of the FI phase in the Pr/Ca system with x ~ 0.2 while 
the La/Sr system with similar hole concentrations are 
ferromagnetic metals. 

In this work, we study a charge and orbitally ordered 



state which can be viewed as an orbital polaron lattice 
and which produces the fully polarized FI state for x 
= 1/4. Hartree-Fock (HF) calculations on d-p-type lat- 
tice models show that the magnetic coupling between the 
orbital polarons is actually ferromagnetic and support 
this idea. We also shortly discuss the evolution from 
the FI state at x = 1/4 to the AFI state at x = 1/2 in 
Pri_2;Ca x Mn03 as well as the relationship between the 
orbital polarons and the small ferromagnetic polarons ob- 
served in Lai_ x Ca x Mn03. 

II. ORBITAL POLARON 

First, let us explain the basic idea pader lying the con- 
cept of orbital polarons (see also refli-3). Since the elec- 
tronic configuration of Mn 3+ in i?Mn03 is t\ g e g , m an 
octahedral coordination, there exists double orbital de- 
generacy and a strong Jahn- Teller effect. On the other 
hand, the e g orbitals are empty in Mn 4+ . When one 
dopes i?Mn03 with holes, namely, puts Mn 4+ ions in 
the background of the Mn 3+ ions, the e g orbitals of all 
the Mn 3+ site surrounding the Mn 4+ site tend to be di- 
rected towards it as displayed in Fig. [I]. Such an orbital 
orientation occurs for two reasons. One is simply steric: 
oxygen ions sitting in between the Mn ions move towards 
the Mn 4+ site and, consequently, the MnOg octahedra of 
the neighboring Mn 3+ sites are elongated along the axis 
pointing to the Mn 4+ site. Another factor is that the or- 
bital occupation helps to optimize the covalency between 
the Mn 3+ and Mn 4+ sites: the e g orbitals directed to- 
wards Mn 4+ site allow for maximal Mn-O-Mn hopping. 
An important consequence of such an orbital ordering 
is that, according to the Goodenough-Kanamori rules, 
the exchange interaction between the Mn 3+ and Mn + 
in this cluster is ferromagnetic. Thus one can treat this 
object simultaneously as a lattice (both breathing-type 
and Jahn- Teller-type) polaron, a ferromagnetic polaron, 
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and an orbital polaron. 

For x = 1/4, one can expect a very natural type of 
ordering of these polarons shown in Fig. |2[ In this po- 
laron lattice, the orbital polarons form a body-centered 
cubic lattice which have two such polarons per unit cell. 
Therefore, one can divide the orbital polaron lattice into 
two sublattices which are simple cubic lattices labeled as 
A and B. In Fig. |[ the shaded (open) circles and or- 
bitals indicate the Mn 4+ and Mn 3+ sites in sublattice A 
(B). In each sublattice, since the neighboring orbital po- 
larons share the Mn 3+ sites, the coupling between them 
is very strong and ferromagnetic. On the other hand, 
it is not trivial whether the coupling between the neigh- 
boring orbital polarons belonging to the different sublat- 
tices is ferromagnetic or antiferromagnetic. In order to 
see whether the magnetic coupling between sublattices 
A and B is ferromagnetic or not and to check the stabil- 
ity of this ordered structure (which is actually a form of 
charge ordering (CO)), we performed HF calculation on 
the perovskite-type lattice model with the Mn 3d and O 
2p orbitals under the lattice distortion shown in Fig. [j]. 

III. MODEL HARTREE-FOCK CALCULATION 

We employ the multi-band d-p model with 16 Mn and 
48 oxygen sites in which full degeneracy of Mn 3d orbitals 
and the O 2p orbitals are taken into account. In this 
model, the intra-atomic Coulomb interaction between the 
3d electrons is considered in terms of Kanamori param- 
eters u, u', j and j'. The charge-transfer energy A is 
defined by e° — e p + nil, where e° d and e p are the energies 
of the bare 3d and 2p orbitals and U (= u — 20/9 j) is the 
multiplet-averaged d—d Coulomb interaction. The trans- 
fer integrals between Mn 3d and O 2p orbitals are given 
in terms of Slater-Koster parameters (pda) and (pdn) 
and those between the O 2p orbitals are expressed by 
(ppa) and (ppn). Here, the ratio (pda) / '(pdir) is -2.16. 
A, U, and (pda) for i?Mn0 3 are 4.0, 5.5, and -1.8 eV, 
respectively, which are deduced from the photoemission. 
studiesEj and ab-initio band-structure calculations J_3't£l 
(ppa) and (ppir) are fixed at -0.60 and 0.15 for the undis- 
torted lattice, which are close to the salues widely used 
for various 3c? transition-metal oxides O 

The tig-tig antiferromagnetic coupling is correctly cap- 
tured by the present method which can reproduce the 
A- type AFI .state for LaMn0 3 and the G-type AFI state 
for LaCr03.Ej When the lattice is distorted, the transfer 
integrals are scaled using Harrison's prescription. In this 
model calculation, A is the most important parameter 
and the other parameters are not sensitive to the calcu- 
lated result. Since A typically has an error bar of ± 1 eV, 
we have performed calculations for A = 2.0 and 6.0 eV 
and confirmed that our conclusion is not changed. In this 
work, the magnitude of the distortion is given by the ratio 
(di —d s )/d where d s is the Mn-0 bond length at the Mn 4+ 
site and di is the longest Mn-0 bond at the Mn 3+ site. 



d = (di +d s ) j 2 is the Mn-0 bond length before the lattice 
distortion is included. Please note that (dda), (ddn), and 
(ddS) between the Mn 3d orbitals have been neglected. 
This is because, in the perovskite structure, the short- 
est distance between the Mn ions is approximately 2d 
(~ 4 A) and that, using the Harrison's relation,EZl (dda), 
(ddk), and (ddS) are ~ -0.08, 0.04, and -0.01 eV which 
are very small compared to (pda) and (pdit). Therefore, 
the effective d band width mainly arises from (pda) and 
(pdn) in the perovskites. 

In Fig. |, the energy difference between the ferromag- 
netic state and the AFI state, where the two sublattices 
are antiferromagnetically coupled, is plotted for x = 1/4 
as a function of the distortion of the oxygen octahedron 
(di — d s )/d. The postulated distortion compatible with 
the orbital polaron is justified by the fact that the AFI 
state with the orbital polaron exists as a meta-stable so- 
lution even without the lattice distortion. Without the 
lattice distortion, the FM state is lower in energy than 
the AFI state. The distortion larger than 0.075 opens 
a band gap for the ferromagnetic state and induces the 
FM to FI transition. As plotted in Fig. ||, the magnitude 
of the band gap increases monotonically with the lattice 
distortion exceeding the critical one ~ 0.065. The mag- 
nitude of the charge ordering, namely, the difference of 
the occupation between the Mn 3+ and Mn 4+ sites (ANd) 
also increases with the lattice distortion (AJVj ~ 0.7 for 
(di-d s )/d = 0.1). 

The FI state found here is exactly the orbital polaron 
lattice shown in Fig. |[ Interestingly, under the lattice 
distortion, the FI state is still lower in energy than the 
AFI state. This indicates that the magnetic coupling be- 
tween two orbital polarons in different sublattices is also 
ferromagnetic. In the FI state, the oxygen between the 
two Mn 3+ sites has hole concentration of ~ 0.2 which 
is as large as that in the oxygen between the Mn 3+ and 
Mn 4+ sites. Namely, the holes at the Mn 4+ site are par- 
tially transferred to the oxygen between the two Mn 3+ 
sites and make the Mn 3+ - Mn 3+ coupling ferromagnetic. 
Therefore, the Mn 3+ - Mn 3+ coupling at x = 1/4 is dif- 
ferent from that at x = 0. In this case, once the band 
gap opens due to the lattice distortion, the energy dif- 
ference between the FI and AFI states is not sensitive to 
the magnitude of the distortion. Here, it should be noted 
that the FM state without lattice distortion is homoge- 
neous and metallic, indicating that the lattice distortion 
is essential to stabilize the orbital polarons and realize 
the CO FI state. The orbital polaron lattice, namely, 
the FI state with the polaron-type charge and orbital or- 
dering might be responsible for thc—FI phase found in 
Pri^Ca^MnOa for 0.1 < x < 0.31TB In the La/Sr sys- 
tem, the lattice distortion is suppressed and the orbital 
polaron is not formed. On the other hand, in the Pr/Ca 
system, the lattice distortion is favored and the orbital 
polaron lattice would be stabilized. 
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IV. COMPARISON BETWEEN CHARGE 
ORDERINGS AT X = 1/8, 1/4 AND 1/2 

The charge and orbital ordering in the icy-plane of the 
orbital polaron lattice is shown in Fig. [|(a) for x = 1/4. 
One can visualize it as the orbital zigzags constructed 
from the 3x 2 — r 2 /3y 2 — r 2 orbitals and the Mn 4+ sites 
which are arranged in a staggered way, so that the neigh- 
boring zigzags (shown by the thick solid and broken lines 
) meet at the Mruj + sites. In the C-E-type AFI state 
for 0.3 < x < 0.5,0 the non-crossing zig-zag is_s,tabilized 
by the Jahn- Teller distortion at the Mn 3+ sitetj and the 
charge and orbital ordering along the zigzagstil is shown 
in Fig. ||(b). Thus this orbital ordering is quite differ- 
ent from the one suggested above for x — 1/4. It would 
be interesting to see how the orbital polaron lattice for 
x = 1/4 evolves to the C-E-type AFI state with charge 
and orbital ordering for x = 1/2 which consists of non- 
crossing orbital zigzag: via continuous flipping of "bro- 
ken" zigzags of Fig. H(a) or as a first order phase transi- 
tion with eventual phase separation into "1/4" -like and 
"1/2" -like phases. 

It is instructive to describe the difference between 
two FI states which the model HF calculations give foe 
Pr 3 / 4 Cai/4Mn03 (this work) and for Lar/gSri/gMnOs.Efl 
The calculations for x = 1/8 show that some ferro- 
magnetic states with orbital and charge modulations 
are stable even without a lattice distortion. The pre- 
dicted FI state with charge and orbital ordering for 
Lar/gSrjygMnOs has the hole-rich and hole-poor planes 
alternating along the z-axis. Since the orbital ordering 
in the hole-poor plane is essentially the same as that in 
LaMnOa, a .small amount of lattice distortion can open 
a band gap.E3 The calculated results qualitatively agree 
with the experimental facts that Lai-^Sr^MnOa of x = 
1/8 is a FI state with orbital,-or(dering but does not have 
substantial lattice distortion.ExlLj On the other hand, at 
x = 1/4, the homogeneous FM state is very stable and 
no FI state with orbital and charge modulations is ob- 
tained unless a large lattice distortion is included. In 
the FI state obtained for P^^Ca! /4Mn03, the large lat- 
tice distortions of the Jahn- Teller and breathing type are 
required to open a band gap. Pr 3 /4Cai/4Mn03 can be 
insulating probably because the large-buckling of the Mn- 
O-Mn bonds ( Z Mn-O-Mn ~ 155°H) reduces the band 
width of the e g band and allows the oxygen ions to relax 
easily. In this sense, the FI state in Lay/gSrjygMnOs is in 
the weak coupling limit and is different from the FI state 
for P^^Ca^MnOa which requires the strong electron- 
lattice coupling. La3/4Sr!/4Mn03 however remains a fer- 
romagnetic metal because the straight Mn-O-Mn bonds 
make the bandwidth relatively large and increase the fre- 
quency of the Mn-0 stretching mode. 

It is expected that the electron-lattice coupling in- 
creases in going from the La/Sr to La/Ca to Pr/Ca sys- 
tems. This is consistent with the fact that the C-E-type 
AFI state becomes stable in going from La^Srj / 2 Mn03 



to La^Ca^MnOa to Pr^Ca^MnOa. Therefore, it is 
reasonable to speculate that the orbital polarons still ex- 
ist in metallic Laa^Ca^MnOa but cannot form a lattice 
to realize the FI phase because the electron-lattice cou- 
pling in La3/4Cai/4Mn03 is not strong enough compared 
to Pr 3 /4Ca! /4Mn03. Actually, small ferromagnetic po- 
larons are observed in Lai-xCa^MnOsQ and their size 
is comparable to that of the orbital polaron shown in 
Fig. [l| The tendency to form orbital polarons. discussed 
above may play a role in the phase separation.!!] In partic- 
ular, the important role of lattice distortion which follows 
from our results (see Fig. [|) may in principle help to ex- 
plain the large scale phase separation observed in ref.B: 
In different degree of lattice distortion in different parts 
of the sample may stabilize respectively the homogeneous 
metallic phase in one of them and the CO insulating state 
in another, the charge densities in them being (almost) 
equal. 

At the end of this section, let us discuss possible limita- 
tion of the present calculation. First, we cannot exclude 
the possibility that other lattice distortions would sta- 
bilize some charge and orbital orderings which are not 
considered here. For example, the Jahn- Teller distortion 
of LaMn03-type would stabilize the charge-and orbital 
ordering similar to that predicted for x=l/StB. With the 
hole-poor plane having the orbital ordering of LaMn03, 
the hole-rich plane would have the hole concentration x 
— 1/2 and may have the CE-type state. However, it is 
not clear whether such layered structure would support 
a ferromagnetic state or not. Second, the Pr 4f orbitals 
are not included in the present model. It is possible that 
the Pr 4/-0 2p^ band accommodates some holes and, 
consequently, the effective hole concentration in the Mn 
3d-0 2po- band is reduced. In order to check this possi- 
bility, one should study a model with the Pr 4/ orbitals 
in future. 



V. CONCLUSION 

In conclusion, the model HF calculation indicates that 
orbital polarons, in which the orbitals of Mn 3+ sites are 
directed towards the Mn 4+ -like site and the spins are fer- 
romagnetically aligned, are stabilized by the breathing- 
type and Jahn- Teller-type lattice distortions and that the 
magnetic coupling between these polarons is ferromag- 
netic. At x = 1/4, they form a body-centered cubic lat- 
tice (a form of charge ordering), which can be also viewed 
as a staggered arrangement of the orbital zigzags. Fur- 
ther theoretical and experimental studies are required in 
order to check whether the proposed structure is realized 
in Pr 3 /4Cai/4Mn03, and in general to clarify the impor- 
tance of the orbital polarons and their relation to the 
CMR effect and to the phase separation. 
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FIG. 1. Schematic drawings of the orbital polaron around 
Mn 4+ -ion. The shaded orbitals and circles indicate the Mn 3+ 
and Mn 4+ sites, respectively. The open circles show the oxy- 
gen ions sitting between the Mn 3+ and Mn 4+ sites. The ar- 
rows indicate the shifts of the oxygen ions. 




FIG. 2. Schematic drawing of the orbital polaron lattice 
for Pr 3 /4Ca 1( / 4 Mn03. The shaded (open) orbitals and circles 
indicate Mn 3+ and Mn 4+ sites in sublattice A (B). 
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FIG. 3. Energy per formula unit cell of the ferromagnetic 
state (Ef) relative to the antiferromagnetic and insulating 
state (Eaf) as a function of the lattice distortion. The mag- 
nitude of the band gap for the ferromagnetic state is also 
plotted as a function of the lattice distortion. In the ferro- 
magnetic (antiferromagnetic) state, spins in sublattice A are 
parallel (antiparallel) to those in sublattice B. 




(a) polaron-type (FI) 
staggered zigzags 




(b) non-crossing zigzags (CE-typeAFI) x=l/2 
FIG. 4. Charge and orbital orderings in xy-pl&ne (a) for 
P r 3/4Ca 1/ / 4 Mn03 and (b) for Pr^Ca^MnOs. Shaded cir- 
cles indicate Mn 4+ sites and open circles indicate Mn 3+ sites 
with the 3z 2 — r 2 -type orbital. In Pr3/ 4 Cai/4Mn03, the neigh- 
boring zigzags shown by solid and broken lines meet at the 
Mn 4+ site. In Pr^Ca^MnOa, the zigzags shown by the 
solid lines never meet each other. 
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